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We re-analyze data from available finite-temperature QCD simulations near the chiral transition, 
with the help of Chiral Random Matrix Theory (chRMT) . Statistical properties of the lowest-lying 
eigenvalues of the staggered Dirac operator for SU(3) lattice gauge theory with dynamical fermions 
are examined. We consider temperatures below, near, and above the critical temperature T c for the 
chiral phase transition. Below and above T c the statistics are in agreement with the exact analytical 
predictions in the microscopic scaling regime. Above T c we observe a gap in the spectral density and 
a distribution compatible with the Airy distribution. Near T c the eigenvalue correlations appear 
inconsistent with chRMT. 



I. INTRODUCTION 

It is inherent to the lattice approach to quantum field 
theory, that one has to extrapolate from finite lattices, fi- 
nite statistics and non-critical coupling parameters to in- 
finite lattices, infinite statistics and critical points. Since 
the result supposedly is a non-trivial, non-perturbatively 
defined quantum field theory, this process is plagued by 
uncertainties. A typical example of such a situation, 
where all these aspects combine, is the study of the ther- 
mal transition in QCD for small quark masses. One is 
interested in the continuum limit (gauge coupling g — > 0) , 
small or vanishing fermion masses (to — > 0), close to crit- 
ical temperature (T — > T c ) in the thermodynamic limit 
(L x — > oo) — a formidable problem. 

The extrapolations are always based on assumptions 
on the asymptotic behavior. Well known examples are 
scaling functions based on renormalization group and chi- 
ral perturbation theory — an expansion around a ground 
state with Goldstone bosons. Here we will examine an- 
other such approach, which should allow the extrapola- 
tion to infinite volume and vanishing fermion mass: Chi- 
ral Random Matrix Theory (chRMT). 



A. Chiral Random Matrix Theory 

RMT attempts to identify universal features of ensem- 
bles of (random) matrices with common symmetry prop- 
erties. Its chiral version, if successful, allows to separate 



two aspects of a theory like QCD: the general universal 
properties shared with other theories from the model- 
specific "dynamical" content of the theory. Microscopic 
eigenvalue distribution shapes are an example for the first 
aspect, expectation values of the fermion condensate for 
the second. 

The limitations for validity of the chRMT considera- 
tions (for a given T in the phase of broken chiral symme- 
try) are set by |],(| 
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< L x < , 

Aqcd mi? 



(1.1) 



where L x is the linear size of the system and is the 
mass of the lightest (pseudo-) Goldstone boson. The sec- 
ond restriction imposes that the pion does not fit into the 
space-time volume and it therefore appears to be unphys- 
ical. However, various correlators in the Dirac operator 
spectrum can be computed precisely in this limit. 

chRMT has been proven to give exact analytical pre- 
dictions for the spectrum of the Dirac operator in the 
microscopic limit ||. The microscopic scaling region is 
simply a blowup of the origin in the spectrum. To be spe- 
cific, one considers eigenvalues A on the scale where 
£ is the chiral condensate, related to the spectral den- 
sity per unit volume p(X) via the Banks-Casher Q re- 
lation, £ = liniA->o hmy^oo 7rp(A). This regime is, by 
definition, only well defined in the spontaneously broken 
phase where p(0) ^ 0. In the phase with restored sym- 
metry the scale of interest is set by the density of states 
in the vicinity of the onset of p(X). 
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Here we present a study of the microscopic correlators 
in the spectrum of the staggered Dirac operator in SU(3) 
gauge theory with dynamical fermions at finite tempera- 
ture. Specifically we examine the low lying eigenvalue 
statistics at temperatures below, near, and above the 
critical temperature of the chiral phase transition. 

Our analysis is based on the evaluation of the MILC 
collaboration's gauge configurations |]||. We therefore 
concentrate on the new aspects connected to RMT ideas 
for the spectral correlators of the Dirac operator. In par- 
ticular for T > T c we study the singularity at the inner 
endpoint of the spectral density. 



B. Temperature transition 



Strictly speaking, in the continuum limit a — > 0, non- 
zero temperature is realized for lattices n x x n t with 



T 



l/(n t a(/3 fl )) and n t /n x 



0. In that limit, for 



vanishing quark mass to, one expects a phase transition 
at T c . In H the critical temperature was estimated to lie 
between 143 and 154 MeV. For T < T c chiral symmetry 
is broken spontaneously, with massless pseudoscalars and 
(tpip) 7^ 0; above T c we expect restoration of this sym- 
metry. Whereas for pure Yang-Mills theory the decon- 
finement transition is associated with a breaking of the 
center-symmetry with the Polyakov loop as order param- 
eter, this symmetry is explicitly broken by the fermion ac- 
tion. Nevertheless, remnants of the original breaking fea- 
ture of the Polyakov loop persist even for small fermion 
masses. 

The nature of the chiral phase transition depends on 
the number of flavors Nf. An argument based on a 3D a- 
models analysis predicts a first order phase transition 
for Nf > 3. For Nt = 2 one expects a second order phase 
transition with SU(2) x SU(2) ~ 0(4) scaling behav- 
ior. However, even first order behavior may be arguable 
|^| Eo[. For staggered fermions at non- vanishing lattice 
spacing the correct counting of flavors is unclear since 
flavor symmetry is restored only in the continuum limit. 
Staggered fermions (as simulated by MILC with the hy- 
brid R-algorithm Q]) correspond to the case Nf = 2 
in the continuum limit. On coarse lattices the so de- 
fined fermions should show at least [7(1 ) ~ 0(2) scaling 
behavior. For a discussion of the various scenarios cf. 
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It is unclear whether the phase transition at T c extends 
towards to > or whether, when moving from lower to 
higher temperature, one just observes crossover-like be- 
havior. Some evidence points towards this second sce- 
nario |T^,|l2|. In the following, we denote the crossover 
(phase transition) position by T c (m) or simply T c . 



II. EXPECTATIONS FROM CHRMT 

According to the nature of the elements in the random 
matrix, chRMT appears in three universality classes. 
In this paper we consider the SU(3) gauge theory with 
quarks in the fundamental representation, which belongs 
to the universality class of the chiral unitary ensemble 
(chUE). The partition function under study is P, |l^ , |l7[ 

Z lNf \{m}) = J dM e - NTrU ^ 

N f 

x Y[ det(M + T + im f ) , (2.1) 
/=i 

where M is a 2N x 2N block hermitian matrix (the ele- 
ments of W being random complex numbers), and T is a 
deterministic, i.e. not-random, off-diagonal block matrix 



M 



W 




T 



t 
t 



(2.2) 



Here dM denotes the Haar measure, U(x) is an analytic 
function. 

The predictions from chRMT concern the correlations 
between the eigenvalues A of D = M + T on the scale 
of individual eigenvalues in the thermodynamic limit 
N — > oo. The matrix D is the analogue of the mass- 
less Dirac operator in QCD. The chiral phase transition 
within chRMT is identified through the value of the spec- 
tral density of the eigenvalues of D near zero, i.e. using 
the Banks-Casher relation. 

Modeling the chiral phase transition in chRMT 
amounts to driving a depletion of eigenvalues of D near 
the origin by means of some temperature parameter. 
Two separate approaches have been examined in the lit- 
erature. First , th e unitary invariant chRMT ^|J|], corre- 
sponding to ( |2.1[ ) with T = 0, in which the chiral phase 
transition is driven by tuning U(M 2 ). Second, the non- 
unitary invariant chRMT 16-|l8|, corresponding to (2.1) 
with U(M 2 ) = M 2 , where the deterministic block matrix 
T mimics the effect of the temperature. In this paper we 
do not need to distinguish between the two approaches 
as they are consistent for the quantities measured here. 

Below T c , i.e. when p(0) ^ 0, chRMT predicts ||J the 
probability distribution for the smallest eigenvalue (for 
the trivial topological sector) 



P 



, r = £ e - z 2 /4 deti^-g^ C, 3 ({Vm 2 +z 2 }) 
2' &eti<i,j<N f Mi ({a*}) 



(2.3) 

with 

C»j(M) = Mr lj i+i(^) , 

z = 2ttX P (0)N , and fi f = 2irm f p(0)N . (2.4) 
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where p(0) is the spectral density at the origin for the 
massless situation (i.e. when m/ = in ( |2.l|) ). Ij de- 
notes the j'th modified Bessel function. This result is 
universal in the chRMT context, that is, the analytic 
form of P(z, {p}) does not change under deformations of 
U(M 2 ) provided that p(0) ^ 0. After the identification 
V = 2 N (V is the physical volume in lattice units), ( |2.3| ) 
allows to extract £ = tt p(0) (the fermion condensate in 
the chiral limit) from finite-volume Dirac spectra. Of 
course, in this the mild condition (1.1) must be satisfied. 

Above T c , when there is a finite gap in the spectral 
density, the repulsion between the eigenvalue pair ±A m i n 
becomes negligible; chRMT hence predicts a soft inner 
edge, known as the Airy-solution [ pOpl] ] . 

At T c — signaled in chRMT by a power-like behavior 
of the spectral density at small A — p(X) oc A« , the pre- 
diction from chRMT is not unique. It turns out that the 
spectral correlators depend on the value of 6 p2] ]. 

The distribution of the smallest eigenvalue is a spec- 
tral one-point correlation function and is quite sensitive 
to statistical fluctuations (see below). As an additional 
measure we also study a two-point correlator: the level 
spacing distribution P(s). Note that the level spacings, 
s = Si+i — Si, are determined in the unfolded spectrum 
{si}i=i- Unfolding separates the fluctuation properties of 
the spectrum from the supposedly smooth background 
behavior. The unfolded variable is defined in terms of 
the eigenvalue spectrum and the local average spectral 
density by 



(p(A))dA 



(2.5) 



The RMT prediction for the level spacing distribution is 
well approximated by the unitary Wigner surmise 



(2.6) 



The level spacing distribution is not expected to be af- 
fected by temperature and masses in chRMT, see e.g. 
II- 

chRMT makes predictions for average spectral corre- 
lators in sectors with definite topological charge v, i.e. 
derived assuming exact zero modes (these are not in- 
cluded in the predicted distributions; cf. for the re- 
sult of the weighted summation of all topological sec- 
tors). For Ginsparg- Wilson fermions [B5[ , which realize 
chiral symmetry on the lattice, one may identify exact 
zero modes as resulting from topological excitations ac- 
cording to the Atiyah-Singer index theorem (for Wilson 
fermions one can hypothesize that zero modes are re- 
placed by real modes). This is not the case for staggered 
fermions, where exact zero modes are absent away from 
the continuum limit p6| and even gauge configurations 
with non-vanishing topological charge do not give zero 
eigenvalues. Exact zero modes are here replaced by "al- 
most" zero modes which accumulate to the origin in the 



continuum limit. In the strong coupling region the mi- 
croscopic staggered Dirac spectra summed over all topo- 
logical sectors show [^,|8| good agreement with the an- 
alytical prediction for the topologically neutral, v = 0, 
sector from chRMT. However, approaching weaker cou- 
pling observations contradicting this scenario have been 
found in a two-dimensional context p9| . 

Before turning to the numerical studies, let us com- 
ment on the validity of the chRMT predictions. The 
condition for application of chRMT in lattice analyses is 
well established when T ~ 0: The range in the unfolded 
spectrum over which the chRMT correlations dominate is 



|A|< 



nl, where is the pion decay constant. 



An equivalent statement for T > T c is not known and no 
stringent tests of the low-lying eigenvalue statistics have 
been carried out so far. 

Let us emphasize that even though the larger part of 
the studies of chRMT have been focused on the situa- 
tion where p(0) ^ 0, there is nothing wrong from first 
principles in using chRMT when p(0) = 0. 



III. GAUGE FIELD CONFIGURATIONS AND 
ANALYSIS 

By courtesy of the MILC collaboration Q there are 
sets of gauge configurations |(| available to the lattice 
community. These were generated with two species of dy- 
namical staggered fermions, at various lattice sizes, tem- 
peratures, values of the gauge coupling and small values 
of the bare fermion mass. In Table | the samples used in 
the present study are listed. For further details on the 
method of determination of the gauge configurations and 
the physical parameters we refer to pi. 



TABLE I. Summary of the MILC configurations used in 
our analysis (for FT01 we only considered a subset of the total 
of 149 configurations available). The suggestions in the right- 
most column are based on MILC's results. The transition is 
near /3=5.26 for n t = 4 and /3=5.725 for n t = 12. 



MILC-set 


# conf's 


n x 


n t 


P B 


ma(p g ) 


phase 


124 A 


61 


12 


4 


5.25 


0.0125 


< T c 


124B 


91 


12 


4 


5.26 


0.0125 


< T c 


124C 


126 


12 


4 


5.27 


0.0125 


near T c 


124D 


57 


12 


4 


5.28 


0.0125 


> T c 


1241 


42 


12 


4 


5.25 


0.008 


< T c 


1242 


50 


12 


4 


5.255 


0.008 


< T c 


1243 


45 


12 


4 


5.26 


0.008 


near T c 


1244 


47 


12 


4 


5.265 


0.008 


near T c 


1245 


42 


12 


4 


5.27 


0.008 


> T c 


1246 


40 


12 


4 


5.28 


0.008 


> T c 


FT01 


30 


24 


12 


5.65 


0.008 


< T c 


FT03 


131 


24 


12 


5.725 


0.008 


near T c 


FT04 


188 


24 


12 


5.8 


0.008 


> T c 


FT05 


146 


24 


12 


5.85 


0.008 


> T c 
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The massless staggered Dirac operator is anti- 
hermitian and therefore has purely imaginary eigenval- 
ues, lying symmetric to the origin. We determine the low- 
est lying eigenvalues with help of the implicitly restarted 
Arnoldi method ]3^ | , using Chebyshev polynomials to im- 
prove the convergence. 

The convergence speed depends on the separation bet- 
ween the eigenvalues. The configurations below T c there- 
fore exhibit much slower convergence and for V = 24 3 x 
12 the determination is then quite time-consuming. Con- 
figurations that are supposedly above T c develop a gap 
for the smallest eigenvalue and the convergence of the 
diagonalization is faster. All eigenvalues have been ob- 
tained with a precision of at least 5 significant digits. 

We have determined (on the positive imaginary axis) 
the lowest ten eigenvalues for the 12 3 x 4 ensembles and 
the lowest eight eigenvalues for the 24 3 x 12 ensembles. 

The chRMT-prediction for the smallest eigenvalue dis- 
tribution has been tested against quenched QCD | f33|]34| 
and dynamical SU(2) |2£| lattice simulations at T = 0; 
there are also recent quenched QCD results for non- 
zero temperature |3^]. RMT distributions for the bulk 
nearest-neighbor spacing in both, the confinement and 
the deconfinement phase of full QCD have been obs erv ed 
in p3] on 6 3 x 4 lattices (Wigner surmise, see Eq. (2.6)). 



IV. REPORTING THE RESULTS 

We now turn our attention to the measurement of the 
statistical properties of the Dirac operator. Starting with 
the distribution of all eigenvalues and then separating out 
the lowest and the second lowest eigenvalue we observe 
the first indication of a discrepancy with chRMT. In or- 
der to investigate this discrepancy further we study the 
level spacing distribution. Finally we focus on the in- 
ner edge of the spectrum for T > T c and measure the 
exponent for the singularity of the spectral density. 



A. Eigenvalue density 

Fig. [l] gives the eigenvalue distribution as obtained 
from the lowest 10 (or 8 for the large lattices) eigenvalues 
for each configuration. It coincides with the distribution 
from all eigenvalues only up to the lowest of all 10th (8th) 
eigenvalues. In order to emphasize this feature we also 
plot the distribution histogram of this 10th (8th) eigen- 
value in full black. 

A remarkable change in the features of the distribu- 
tions occurs around T c , as given in Table |. In particular 
for n x = 24, nt = 12 a gap in p(X) appears for T > T c . 

Since we do not know about the topological charge of 
most of the configurations studied here, we cannot sepa- 
rate the trivial from the non-trivial topological sectors, as 
would be required for a faithful comparison with chRMT 
distributions. For sufficiently rough lattices (i.e. in the 
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FIG. 1. Histograms for the 10 (or 8) smallest eigenvalues. 
The contribution from the 10th (8th) eigenvalue is indicated 
in black, (a) 12 3 x 4, m = 0.0125, (b) 12 3 x 4,m = 0.008, (c) 
24 3 x 12, m = 0.008. 



strong coupling region) , one can argue |2S| that the topo- 
logical charge v is effectively zero from the fermionic 
point of view; however this is no more the case for fine 
enough lattices |^| and the problem of the knowledge of 
the topological charge becomes critical. For T > T c the 
situation is completely different, since topological fluctu- 
ations are suppressed in the continuum theory. 
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FIG. 2. Histograms for the 1st eigenvalue. (a) 
12 3 x 4,m = 0.0125, (b) 12 3 x 4,m = 0.008, (c) 
24 3 x 12, m = 0.008. For the data where we fitted to 
chRMT-distribution (cf. Table II) we also plot the fits. The 
error bars in (a) are shown to indicate the typical size of errors 
in all histograms. 



The "microscopic" distributions should be in the 



chUE universality class and prediction (2.3) should ap- 
ply in particular for the smallest eigenvalue for T < T c . 
A fit of the corresponding prediction for the topologically 
trivial case (y = 0) to the normalized data provides us 
with the (infinite volume) parameter S, We assume the 
continuum symmetry, i.e. Nt = 2 in (2.3). This one- 
parameter fit appears reasonable only for data concern- 
ing temperatures well below T c as indicated in Fig. |](a) ■ 



This is made explicit in Table || where the fitted values 
of £ together with the corresponding ^/d.o.f. are re- 
ported; the latter increases with T, an d for T ~ T c (in 
agreement with Table I) prediction (2.3) becomes incom- 
patible with data. 



The formal chRMT expression (2.3) gives the eigen- 
value distribution as a function of volume and fcrmion 
mass; the parameter S (the spectral density at the ori- 
gin) is thus defined implicitly as the extrapolation to in- 
finite volume and vanishing fermion mass. From our fits 
- if the data follows chRMT formulas - S should there- 
fore be independent of the spatial lattice size L x and the 
fermion mass m. 

For e.g. 12 3 X 4, p = 5.25 we find agreement of £ for 
rn = 0.0125 and 0.008 within the errors. At P = 5.26 
the values disagree. This value of P is, however, close to 
T c and the position of the phase transition (or crossover) 
is quite sensitive to m. Such a change of the transition 
point is not accounted for in chRMT, which is not at all 
sensitive to the underlying dynamics of QCD. 

For comparison we have to extrapolate the MILC val- 
ues H for the chiral condensate both, to infinite volume 
and to vanishing quark mass. Since (except for chRMT) 
we have no firm prediction concerning the volume depen- 
dence, we just extrapolate the MILC results for the large 
lattices at p = 5.65 and 5.725 linearly to vanishing quark 
mass. We know however that finite- volume effects on the 
condensate increase as the quark mass decreases. There- 
fore, it is not surprising that our linear extrapolation of 
MILC results tends to come out slightly but systemati- 
cally smaller than the values we find in Table ||. 

The issue of the topological sector is likely to be par- 
ticularly relevant for the finest lattice at our disposal 
(24 3 x 12), where almost zero modes could be present 
and spoil the validity of the trivial sector predictions from 
chRMT. These could be the cause of the bump observed 
for small A at P = 5.8, both in the spectral density and in 
the smallest eigenvalue distribution. It is therefore with 
some hesitation, that we compare the histograms for the 
smallest eigenvalues in Fig. H with these predictions. 



TABLE II. Values of the scale £ as obtained from compar- 
ison of lattice data for P(A m i n ) with chUE. 



n x 


n t 


P 


m 


E 


X 2 /d.o.f. 


12 


4 


5.25 


0.0125 


0.647(29) 


0.495 


12 


4 


5.26 


0.0125 


0.571(26) 


0.562 


12 


4 


5.27 


0.0125 


0.449(29) 


1.651 


12 


4 


5.28 


0.0125 


0.138(13) 


2.129 


12 


4 


5.25 


0.008 


0.686(37) 


0.368 


12 


4 


5.255 


0.008 


0.551(52) 


1.218 


12 


4 


5.26 


0.008 


0.459(43) 


1.379 


12 


4 


5.265 


0.008 


0.331(39) 


2.517 


12 


4 


5.27 


0.008 


0.168(21) 


2.806 


24 


12 


5.65 


0.008 


0.0198(18) 


1.048 


24 


12 


5.725 


0.008 


0.0063(60) 


5.063 



5 



In order to further investigate this feature, we studied 
(for the set with (3 = 5.725) the influence of the config- 
urations where the eigenvectors uq of the lowest eigen- 
values exhibit a large contribution to the total chirality, 
i.e. |(uo75Mo}| > 0.08. According to the index theo- 
rem, these configurations with large chirality, which make 
up roughly one half of the ensemble, tend to carry non- 
vanishing topological charge and therefore zero modes. 
Indeed we find that a substantial part (75%) of the first 
peak may be explained from these contributions. 

These findings suggest that indeed topological modes 
are responsible for a low- lying peak in the distributions. 
Below T c all topological sectors are present and the low- 
lying eigenvalues have comparable magnitudes (their av- 
erage position being roughly proportional to i>). When 
the temperature approaches T c the theoretical expecta- 
tion is that the topological fluctuations begin to be sup- 
pressed, although still present in the ensemble, quasi su- 
perimposed on the background distribution, which starts 
to broaden significantly with increasing temperature. 
Sufficiently far above T c only the topologically trivial 
sector survives and there may be no small eigenvalues 
at all. This is indeed what we actually observe for the 
lattice 24 3 x 12. 

In a recent study of quenched configurations there 
have been indications for a dilute gas of instanton-anti- 
instanton pairs producing a Poissonian distribution of 
small eigenvalues above T c . These may be suppressed 
or absent when considering dynamical fermions. In our 
context this seems to be the case for the finest lattices 
(24 3 x 12 with (3 = 5.85) at our disposal. 

Fig. U gives the histograms for the 2nd smallest eigen- 
values. Again we notice a dramatic change of the distri- 
bution shape around T c (m). 
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FIG. 3. Histograms for the 2nd smallest eigenvalue for 
12 3 x 4, m = 0.008. By comparison with Fig. 0(b) we ob- 
serve that the mutual overlap between P(A m i n ) and P(\2) 
increases as /3 — > (3 C . This is not consistent with chRMT. 

We interpret the sudden flatness of the distribution of 
the smallest eigenvalues as being (i) due to the vanish- 
ing spectral density and (ii) due to increasing statistical 
fluctuations near the chiral phase transition. The latter 



effect is not reproduced in chRMT since this is a zero di- 
mensional and non-dynamical theory. Furthermore, the 
mutual overlap of P(\x) and PiM) increases for (3 — > f3 c . 
This is also inconsistent with chRMT. In order to study 
this effect further we now turn to the level spacing dis- 
tributions. 



B. Level spacing distribution 

Another observable with definite predictions from 
RMT (cf. Sect. ||) is the distribution of level spac- 
ings. The advantage here is, that the level spacing should 
not be influenced by possible distortions of the smallest 
eigenvalues due to the unknown topological charge of the 
configurations (if the smallest eigenvalues are removed 
from the data). 

The studies of the level spacing statistics in lattice data 
so far have shown a u niform picture consistent with the 
RMT prediction (2J3). The agreement extends on both 
sides of the confinemcnt-dcconfinement phase transition 
[^3|,^6|. However, to the authors' knowledge there is no 
analytical prediction from chRMT for the level spacing 
distribution when focusing on the soft edge or at T c . So 
one might worry that the standard prediction, Eq. (|T|), 
is not appropriate when p(0) — 0. Within the T-model of 
|l6| for T = T c and T = 3 T c , we have performed a numer- 
ical high statistics simulation to elimina te s uch doubt, 
and we there confirmed the distribution ( |2.6|) . (The T- 
model is defined by ( |2~i"| ) with U(M 2 ) = M 2 and t in 
(|2.2[) chosen proportional to the unit matrix.) 



Usually it is possible to get high statistics on the level 
spacing distributions since each configuration provides a 
large number of eigenvalue spacings. However, already 
the lowest 10 (or 8) eigenvalues allow for a crude esti- 
mate of the distribution shape. 

Recall that the level spacing dis tribution is measured 
in the unfolded spectrum, see ([D]). Here we use the av- 
erage spacings (Ai + i — Xi) between contiguous eigenvalues 
to define the unfolded level spacings 

(Ai+i - A 4 ) 

s z+ i-s z = — — . (4.1) 

(Ai+i — Ai) 

In Fig. ^ we compare the data with the parameter-free 
theoretical expectation. Whereas below and above T c we 
find reasonable agreement with the theoretical expecta- 
tion, there are clear discrepancies near T c . We observe 
unexpected high histogram entries. Since the average 
value by definition is 1 this then leads to a shift of the 
central peak to the left. 

In order to further check our unfolding procedure, we 
also considere d o ther approaches, e.g. using a average 
density as in (2.5) by smoothing our distribution in var- 
ious ways. We furthermore tried to discard the higher 
lying eigenvalues, e.g. using only the lowest 5 level spac- 
ings or introducing a cutoff near the peak of the distri- 
butions in Fig. ||. In all those checks we found essentially 
the same behavior with discrepancies near T c . 
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FIG. 4. Histograms for the distribution of (a) the first 9 
eigenvalue spacings on the 12 3 x 4 lattice and (b) the first 7 
eigenvalue spacings on the 24 3 x 12 lattice with m — 0.008 
(accordi ng t o [5,6] /3 C ~ 5.7 in this case). The standard pre- 
diction (2.6) for the level spacing distribution is plotted for 



comparison. 

In conclusion of this section, we observe at T ~ T c a 
breakdown of the otherwise universal microscopic spec- 
tral correlations. The dynamics of QCD plays an essen- 
tial role in the phase transition. A RMT model where 
such dynamics is not there fails to account for the in- 
creased fluctuations in the eigenvalue level spacings. 



p(X,u)=K 



A 2 



1 



u+l/2 



0,2,4,.. 



(4.2) 



where if is a known constant. The constant u takes 
the value u = in the generic ch RM T, i. e. without 
fine-tuning the matrix potential in (2.1) (which would be 
necessary in order to obtain higher values of u). This 
corresponds to a square-root-like eigenvalue density near 
A. 

One concern here is to measure u. With the limited 
amount of data available it is not possible to do this based 
on the spectral density only. Instead we propose to study 
the average distance between the smallest and the sequel 
eigenvalues = (Aj+i — Ax); the extraction of u is 

carried out by noting the following scaling relation in the 
index i, ordering by size the eigenvalues which follow the 
smallest, 



(4.3) 



This proportionality follows by integration in Q4. 

In Fig. H we display the seven average distances, 
from the ensemble of 146 configurations on a 24 3 x 12 
lattice for (3 g — 5.85. We also exhibit the best fit to fl4.3| ) 
with respect to u, giving u — 0.117(71). Also shown are 
the corresponding curves for u = and u = 2. The value 
u = is clearly favored. 
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C. Soft inner edge 



FIG. 5. The behavior of the average distance from the 
smallest eigenvalue to the i'th eigenvalue for the 24 3 x 12 
ensemble at /3 — 5.85. 



We now turn to the results for T > T c . In our results 
for 24 3 x 12 at (i = 5.85 (Fig. |l|(c)) a gap in the spectral 
distributions is obvious. However even at = 5.8 we may 
speculate, that a clear signal of a gap is only prevented 
by the (topological) quasi-zero modes responsible for the 
small bump at small eigenvalues. 

Recall, that chRMT @ predicts the presence of a gap 
in the spectral density p(X) of the Dirac operator cen- 
tered around A = 0. Furthermore, the inner edge of this 
gap is predicted to show a singularity, at a point A, in 
the macroscopic spectral density J37| 



Since the two-point correlations behave as expected 
from RMT, we now turn to the one-point distribution. 
The microscopic behavior of the spectral density in the 
vicinity of this singularity is universal in the chRMT 
sense, but depends on the value of u p7[ . For u = 
the exact analytical prediction for the microscopic spec- 
tral density in the vicinity of the inner edge is p0| ] 

PAi(z,0) = (Ai'(^) 2 -z(Ai(z)) 2 . (4.4) 

Here the origin has been moved to the inner spectral end- 
point A by means of the rescaled eigenvalue z, which is 
defined through 
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X = A 



z 

1 + 2 V TrAis: 



(4.5) 



The consistency with the prediction (4.2) for u = 
and the approximate validity of chRMT correlations in 
the level spacing statistics above T c suggest that the Airy- 
density (4.4) corresponding to the value u = should fit 
the spectral density. If it does, then we can extract the 
inner-endpoint A of the spectrum in the thermodynamic 
limit, by fitting pAi(z,0) with respect to A to the lowest 
part of the spectral density, see Fig. 0. 
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FIG. 6. The u = Airy density with endpoint 
A = A — 0.0205. Approximately the domain of the first three 
eigenvalues is shown. 

This fit does not convincingly confirm Airy density. 
However, the statistical fluctuations at this /3-value af- 
fect the one-point distribution substantially and prevent 
a decisive comparison. 



The observed discrepancy from the RMT level spac- 
ing statistics may also be interpreted as an inclination 
towards Poissonian statistics; distribution shapes inter- 
polating between Wigner and Poissonian statistics have 
been suggested by Brody Q . 

As may be seen from the 24 3 x 12 ensembles at (3 = 5.8 
and 5.85, a gap develops the spectral density for T > T c . 
This is consistent with the observed suppression of topo- 
logical fluctuations in the latter ensemble [^9|. For the 
/3 = 5.85 ensemble we have measured the critical expo- 
nent characterizing the steepness of the density at the in- 
ner edge. The value is found to be compatible with 1/2. 
This is exactly as predicted by chRMT where the chiral 
phase transition is manifested by the crossover from the 
Bessel hard edge to the Airy soft edge. The indications of 
the Bessel to Airy scenario are suggestive but simulations 
with extended statistics are needed in order to quantify 
the observation. However, even with low statistics the 
/3-dependence of the distance between e.g. 8th and 1st 
eigenvalue provides an excellent means to identify the 
change of the phase. 

At low j3, on coarse lattices, staggered fermions ap- 
pear to be blind with regard to the topological charge 
of the gauge configurations, and the smallest eigenvalue 
distribution agrees with the chRMT distribution for the 
v = sector. As the lattice becomes finer, topology be- 
comes more relevant. Although this is maybe "good" 
for the continuum limit of staggered fermions, it affects 
unfavorably the agreement with chRMT since the want- 
to-be-zero modes and the non-zero modes have similar 
eigenvalues, and begin to separate only when the non- 
zero modes are pushed to larger values when increasing 
the temperature. 



V. SUMMARY AND CONCLUSIONS 
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